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Abstract. Thanks to the work of V.I. Arnold Sur la geometrie differentielle 
des groupes de Lie de dimension infinie et ses applications a I'hydrodynamique 
des fluides parfaites (1966) we know that an ideal fluid can be understood as a 
Lagrangian mechanical system on a Lie group of volume preserving diffeomor- 
phisms. Moreover, one can reduce by the particle relabeling symmetry of the 
system to arrive at equations of motion on the Lie algebra of divergence free 
vector fields. In this article we will investigate the equations of motion for a 
flexible body immersed in an incompressible ideal fluid. We will find that the 
configuration manifold for such a system is the source fiber of a Lie groupoid. 
Moreover, we will reduce by the particle relabeling symmetry of the system to 
obtain equations of motion on the Lie algebroid in close analogy with the work 
of V.I. Arnold. As this groupoid acts transitively on itself, we will be able to 
execute this reduction as an instance of Lagrange-Poincare reduction. We will 
also address the case of Navier-Stokes fluids by including a dissipative viscous 
force and amending the original Lie groupoid to include a no-slip condition. 
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1. INTRODUCTION 

Since the work of |Arn66j, it is well known that the dynamics of fluids can 
be described as an ODE on a Lie group. Fluid-structure interactions present an 
interesting challenge because the structure immersed in the fluid interferes with 
the group properties of the system. In this paper we will preserve as much of this 
group structure as possible by describing the system on a Lie groupoid. Compelling 
reasons for doing this are motivated by history. In particular, the publication of 
|Arn66j was soon followed by: 

• proofs of local-time existence of solutions for inviscid fluid equations [EM70J 
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• the creation of computationally tractable Lie group integrators for fluids 
[GMP+11] 

• the geometrization of complex fluids |Hol02[ GBR09J 

• a new paradigm of analysis for the Lagrangian stability of fluid motions 
|Luk88llMis93] . 

See |AK92j and [MR99, §13.7] for an overview. The main theorem of this paper 
parallels that of [Arn66|, and we hope to promote a similar research programme for 
fluid-structure problems. Some specific advantages of having a geometric deriva- 
tion of fluid-structure equations are: 

• We do not make any assumptions on the nature of the immersed body, 
which could be a rigid or flexible medium, or even another fluid. 

• Along the way, we obtain a reduced variational principle for fluid-structure 
interactions. Such principles are a first step in constructing variational- 
geometric numerical integrators (see (HLW021 iMWOTl IMV911 [BS99] ). 

• Having a geometric understanding of fluid-structure interactions might 
prove to be useful in the context of nonlinear stability and control the- 
ory. In |Jacl2j . for instance, the geometric framework established in this 
paper is used to demonstrate the existence of (approximate) limit cycles 
for bio-locomotion. 

In summary, a diversity of perspectives is often beneficial. We hope the content 
of this paper provides a new perspective on fluid structure interaction, and opens 
the door for new types of collaboration between differential geometers, engineers, 
and applied mathematicians. 

1.1. Main contributions. The over-arching contribution of this paper is an un- 
derstanding of fluid-structure mechanical systems using natural operations on Lie 
groupoids. Many of these ideas will look different from the traditional approach 
to fluids. Therefore, this paper also provides a "sanity check", where we make 
sure that the equations of motion in our geometric framework correspond to the 
standard equations of motion in terms of vector-calculus. In particular, the main 
contributions can be listed as: 

• We will identify the configuration manifold for fluid-structure systems as 
the source-fiber of a Lie groupoid. 

• We will identify the particle relabeling symmetry of a fluid-structure sys- 
tems as the isotropy group of the same Lie groupoid. 

• We will state the equations of motion as on ODE on a Lie algebroid. 

• We will verify that the equations of motion stated correspond to the equa- 
tions of motion used in the engineering literature. 

In the following subsections we will review the equations of motion for a structure 
immersed in a Newtonian fluid and provide a preview of the Lie groupoid which 
we will be studying throughout the rest of the paper. 
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1.2. Fluid-structure interactions. The equations that determine the motion 
of a coupled fluid-structure system consist of the Navier-Stokes equation together 
with a forced Euler-Lagrange equation for the evolution of the body. Specifically, 
these equations are given by 

(i) £ (^) - = f. 



Dt\db J db 
du 

(2) — + U -Vu = uV 2 u - Vp, V ■ u = 

L/C 

Here, z/ is the viscosity coefficient, while u and p are the Eulerian velocity and the 
pressure field, Lb is the Lagrangian of the body and fs is the force exerted by the 
fluid on the body, given by the expression for virtual work 

(3) f B -5b = - [ 5b-T-ndS 

JdB 

where T = —pi + v {\Ju + Vw T ) is the stress tensor of the fluid and n{x) is the 
unit external normal to the boundary of the body. For the case of a rigid body 
immersed in a viscous fluid, these equations can be found in, for instance, [MZ05J. 
For the case of general structures immersed in viscous fluids, one can deduce these 
equations from the standard theory of continuum mechanics using the stress tensor 
T, jBatnni §3.3, §5.14]. 

In this paper, we will consider both the case of a viscous fluid interacting with a 
body as well as the inviscid case, for which v — in pj. The boundary conditions 
will be chosen accordingly and are discussed below. 

1.3. Lie groups and groupoids in fluid mechanics. In his foundational work, 
|Arn66j showed that the motion of an inviscid, incompressible fluid moving in a 
fixed container M describes geodesic curves on the group SDiff(M) of volume- 
preserving diffeomorphisms of M. As mentioned before, this geometric insight 
gave rise to a plethora of applications and new results, but is not immediately 
applicable to the case of fluids interacting with other media: as the body moves 
through the fluid, the fluid container changes from one instance to the next, so 
that the fluid configurations can no longer be described by diffeomorphisms of a 
single, fixed manifold. 

A first contribution of this paper is to show that fluid-structure interactions can 
be described by means of a Lie groupoid instead of a Lie group. To explain this 
concept, we recall that the configuration of an immersed body can be described 
by means of embeddings b : B — > lR d where B is a manifold which represents the 
body. The fluid space is the complement of the space occupied by the body, which 
we denote by 

:= closure (R d \b(B)) . 

Now, if the rigid body moves from one configuration, bo, to another, b±, the fluid 
container changes from r^;, to tWb 1 ■ The configuration of the fluid, meanwhile, is 
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given by a (volume-preserving) diffeomorphism ip : p&Jb — > R^y^. This motivates 
us to introduce the set 

(4) G = {(&!,?, &o), &o,6i e Emb(£),<^ G SDiff (r^,^)} 

as the natural configuration space for fluid-structure interactions, where the set 
SDiff (r^6 , PtfUbi) consists of all volume preserving diffeomorphisms from rr^ to 

The set G is equipped with a multiplication, which is only partially defined: two 
elements (6j, (p, b ) and (b^ip', b' ) can be multiplied only if b\ = b' , in which case 

(5) (b[, (f/, b' ) ■ (61, &o) = ¥>' V 2 , M- 

As the multiplication operation is defined only for certain pairs of elements, G is 
not a group like SDiff, but rather a groupoid which is a well studied object within 
differential geometry. For a gentle introduction to Lie groupoids in geometry we 
refer to the survey article [Bro87j . while a comprehensive, geometric overview can 
be found in [Mac05j. More information about Lie groupoids in mechanics can be 
found in [MMdDM06l IWei95j while an application to general relativity using a 
groupoid similar to G can be found in |BFW09| . 




Figure 1. The fluid-structure Lie groupoid G consist of triples 
(&i, ip, bo), where 60, b\ are embeddings from the body reference con- 
figuration into R d , and ip is a diffeomorphism from RRj& to rr^. 
Here is the space available to the fluid when the body is in the 
configuration bo, and similarly for rs^i- 

The infinitesimal version of a Lie groupoid is referred to as a Lie algebroid (see 
Figure [2]), and in this paper, we will show that the Lie algebroid A of Q consists of 
all triples (6, b, u), where b G Emb(i3) and b G T b Emb(£>) are a body configuration 
and velocity, and u is an (Eulerian) fluid velocity field. Both the groupoid G as 
well as its algebroid A therefore have clear and straightforward interpretations in 
terms of fluid dynamics. 

1.4. Derivation of the fluid-structure equations. The second part of our 
paper is devoted to the derivation of the fluid-structure equations ([TJ [2]) using the 
framework of Lagrange-P oincare reduction. In this paper, we start from geodesic 
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Figure 2. Lie groups and pair groupoids are special cases of Lie 
groupoids. Their infinitesimal counterparts are Lie algebras and 
tangent bundles, respectively, and form examples of Lie algebroids. 

motion on the groupoid G and we observe that this system has the group of 
particle relabelings as its symmetry group. This is an extension of Arnold's original 
observation in |Arn66j to the groupoid context. Factoring out by this symmetry, we 
obtain a reduced system, which can be given a particularly convenient description 
once we choose a horizontal subspace of the Lie algebroid. In the last part of the 
paper, we include viscous dissipation into the fluid-structure picture, by modeling 
it as a (non-conservative) force. 

1.5. Previous work. This paper provides a bridge between two areas of study: 
fluid structure interaction and the theory of Lie groupoids. As a result we draw 
from a range of sources on both sides of this divide. In particular we will exploit 
the fact that the Lagrange-Poincare equations can be seen as equations on a Lie 
algebroid. There is much work in the world of fluid structure interaction which 
hints at viewing fluid-structure systems in the light of Lagrange-Poincare reduc- 
tion. Specifically, it is well-known that fluid-structure systems exhibit fluid particle 
relabeling symmetry, and one could surmise that reducing by this symmetry yields 
Lagrange-Poincare equations. However, to the best of our knowledge, this deriva- 
tion does not yet exist in the literature. Nonetheless, similar tasks have been 
executed, and this "gauge theoretic" flavor to fluid mechanics has been adopted in 
a number of instances. For example, |SW89] developed a gauge theory for swim- 
ming at low Reynolds numbers. This gauge theoretic picture of swimming was 
then applied to potential flow in |Kel98] . Building upon |Kel98j . the equations of 
motion for an articulated body in potential flow were derived as Lagrange-Poincare 
equations (with zero momentum) in [KMRMH05J and then expanded to handle 
point vortices in [VKM09J. Similarly, |Rad03j generalized the work of [Kel98] to ar- 
bitrary flows and and adopted the Lagrange-Poincare equations as valid equations 
of motion. 

On the other side of the main theorem for this paper lies the world of Lie 
groupoids and geometric mechanics. The equations of motion we derive will 



ON THE USE OF LIE GROUPOIDS IN FLUID-STRUCTURE INTERACTIONS 



7 



be achieved via a reduction by symmetry of a Lagrangian mechanical system. 
Early work on reduction of Lagrangian systems by symmetry can be found in 
[MS93b, MS93aJ with an intrinsic formulation articulated in |CMR0lj for the case 
of arbitrary finite dimensional configuration manifolds and left-invariant group 
symmetries. Around the same time, a notion of Lagrangian mechanics on Lie 
algebroids was born in [Wei95j and expanded to handle degenerate Lagrangians 
in |Mar01j . It was made clear in |Wei95j that the Lagrange- Poincare equations 
correspond to the Euler-Lagrange equations on a Lie algebroid when the Lie alge- 
broid is transitive. The theory of reduction by symmetries on Lie algebroids was 
developed further in |Mes05] . 

1.6. Organization of the paper. In ^2] we will cover preliminary material to 
be used in proofs and definitions to follow. The reader is encouraged to read this 
section casually and refer back to it as needed. In section §[3] we will define the 
fluid-structure groupoid, G, and its Lie algebroid A. We will also describe the 
configuration manifold Q(b ) as a source fiber of G. Finally we will find that 
A = TQ(bo)/G(bo) where G(b ) is a Lie group. Having laid down the kinematics 
we will be able to construct the Lagrangian on the tangent bundle TQ(bo) in 
section We will find that the Lagrangian exhibits a particle relabling symmetry, 
represented by the Lie group G(b ) and this will suggest that the equations of 
motions can be written on A. In order to perform this reduction from dynamics 
on TQipo) to dynamics on A we will need to understand variations on A better, 
and so we study this in §[5j Using this understanding of variations on A we will be 
able to transfer Hamilton's variational principle and derive reduced equations of 
motion on A in §[6] using Lagrange- Poincare reduction. These equations of motion 
are written in two ways and provide equivalent descriptions of the same system 



from alternative perspectives: one which is more geometric (Theorem 6.1) while 



the other description is the traditional description of fluid structure interactions 



(Theorem 6.2). The proof of Theorem 6.2 is long and relegated to Appendix |B] 
Finally in §[7] we will add viscous dissipation and a no-slip condition to derive the 
equations of motion for a body in a viscous fluid using the same formalism. 
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caring mentor even in death. Many of the ideas in this paper were clarified through 
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the groupoid literature. Additionally, the authors would like to thank Frangois 
Gay-Balmaz, Mathieu Desbrun, Darryl D. Holm, and Tudor S. Ratiu for encour- 
agement and clarifications on fluid mechanics and Lagrange-Poincare reduction. 
Both authors gratefully acknowledge partial support by the European Research 
Council Advanced Grant 267382 FCCA. H. J. is partially supported by NSF grant 
CCF-1011944. J. V. is on leave from a Postdoctoral Fellowship of the Research 
Foundation-Flanders (FWO-Vlaanderen). This work is supported by the irses 
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2. PRELIMINARY MATERIAL 

In this section we will recall some notions from differential geometry which will 
be used in the remainder of the paper. We will assume the reader is familiar with 
manifolds and vector bundles, and we refer to [KN63J IKSM99] for a comprehensive 
overview. All maps will be assumed smooth throughout this paper. Given any 
manifold, M, its tangent bundle will be denoted by tm '■ TM — > M or by TM for 
short. Additionally, given a map / : M — > N, the tangent lift will be given by Tf : 
TM — y TN. We denote the set of fc-forms on M by Q k (M). Given a vector bundle 
tie ■ E -> M we denote the set of ^-valued A;-forms by tt k (M; E) = T(E) ©fi fc (M). 
If a e Q k (M), the exterior derivative of a will be denoted da and similarly we 
may define the exterior derivative of a .E-valued form by d(e <8> a) = e ® da where 
e G T(E). Lastly, pullback and pushforward by a map, tp : M — )■ N, will be 
denoted by (p* and v?* respectively. 



2.1. Vertical bundles. Let 7i"i : E\ — y M and 7T2 : E 2 — >• M be vector bundles 
over M. We define the Whitney sum 

E x © E 2 := {(e x , e 2 ) E Ei x Ezfafa) = n 2 (e 2 )} 

which is itself a vector-bundle over M equipped with the addition operator 

(ei,e 2 ) + (/ 1 ,/ 2 ) = (e 1 + / 2 ,e 2 + / 2 ) 

for any two (ex, e 2 ), (fx, / 2 ) E E x ® E 2 contained in the same fiber. For any vector 
bundle, ti : E — > M, we may invoke the tangent bundle, te '■ TE — > E, and define 
the vertical bundle, ixy : V(E) — > E, where 

V(E) = kernel(TTr) 

and Tiy := te\ V (e)- Equivalently, we may define the vertical bundle as the set of 
tangent vectors of E of the form ^| e _ (e + ef) where /, e are in the same fiber. 
This gives us the following proposition. 



Proposition 2.1. Let v ^ : E © E — > V(E) be the map 

(ex + ee 2 ) . 



v'{ex,e 2 



de 



e=0 



Then is an isomorphism between V(E) and E © E. We call the vertical lift. 
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2.2. Connections. There are three types of connections used in this paper: Levi- 
Civita connections, connections on a vector bundle, and principal connections. 
In this section, we will review the first two notions while principal connections 
will be discussed in Appendix [A] In particular, we define the covariant derivative 
induced by a vector bundle connection, as this concept will be used to formulate 
the equations of motion throughout the paper; see for instance theorems 6A_ and 

Definition 2.1. A horizontal sub-bundle on a vector bundle, n : E — > M, is a 
subbundle, H(E) C TE, such that TE = H(E) © V(E). Given a horizontal sub- 
bundle we define the projections ver : TE — > V(E) and hor : TE — > H(E). The 
projection, ver, is called a connection. 

This definition of connection appears in |KSM99j where it is called a "generalized 
connection" . The choice of a connection has a number of consequences, the most 
important being the existence of a horizontal lift. 

Proposition 2.2. Given a connection on a vector bundle tx : E — >• M , there exists 
an inverse to the map te © Tir\ H , E s : H(E) — > E © TM . We denote this inverse 
by : E® TM -» H(E), and call it the horizontal lift. 

Proof. Choose an e G E and let H e (E),V e (E),T e E denote the fibers of the hor- 
izontal bundle, the vertical bundle, and the tangent bundle above e. As T e E = 
H e (E) © V e (E) and V e (E) = kernel(T e 7r), we see that T e 7r\ H f E s is injective. 
Additionally, T e 7r\ He , E \ must also be surjective since T e n{T e E) = T m M where 
m = 7r(e). Thus T e n\ He{E) is invertible. We define h? : E © TM -> H(E) by 
h^(e,m) = T e Tr\"^j E \ (m). By construction o t e © Tir\ H , E -, is the identity on 

//(/•;). ~ ' " ' □ 

We will occasionally use the notation h\(m) = hJ(e, m). Additionally, the choice 
of a connection induces a covariant derivative. 

Definition 2.2. We define the vertical drop to be the map, v± : V(E) — > E given 
by 

vi = 7r 2 o (v^y 1 

where 7r 2 : E © E — > E is the projection onto the second E component. Let IcR 
be an interval so that we may consider the set C(I,E) consisting of curves from 
I into E. We define the covariant derivative with respect to the connection H(E) 
to be the map, Mi : C(I, E) — > E, given by 



De 



ver 



de 
dt 



Dt 

for a curve e(-) G C(I,E). 

Lastly, the choice of a connection also provides two partial derivatives on a 
vector bundle. 
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Definition 2.3. Let tt : E — > M be a vector bundle and let f : E — )■ R be a 

function. Given a connection we define the partial derivative of f with respect to 
m as the vector bundle morphism, : E — > T*M , given by 

/J^(e),5m\ := (df(e),ht(6m)); 

see |AMR09] . Additionally we define the fiber derivative, ^ : E — > E* , given by 



^(e),e>\ = (df(e),vHe,e>)) := | 



(f(e + ee')). 

e=0 



These partial derivatives allow us to write down the differential of a real-valued 
function in a new way. 

Proposition 2.3. Let f : E -> E. Let be G TE and Sm = Tvr(5e). Then 

W e)>> = (£(e),* m ) + (g<e),^ 

where e(t) is a path in E tangent to be. 

Proof. We observe that 

(df,5e) = (df,hor(5e)) + (df,ver(8e)). 

By definition, the first term is identically (df, h^(T7i(Se))) = {§^, dm}. The second 
term is 

De(t) \ /df De(t) 



(df,ver(5e)) = (df,vl 



Dt ' \de Dt 



□ 



The covariant derivative on curves in E naturally induces a covariant derivative 
on curves in E* given by the condition 

. . d . . / Da \ I De 

(6) Jt M = \-D-v e ! + Y^t 

Remark 2.1. Given vector bundles E 1 and E 2 with horizontal sub-bundles H(E%) 
and H(E 2 ) we can define a horizontal sub-bundle by taking the direct sum: H(E\@ 
E 2 ) = H(Ei) © H(E 2 ). This will be valuable later when dealing with equations of 
motion on the direct sum of a Riemannian manifold and an adjoint bundle. 

2.3. Connections on Riemannian manifolds. The two most important ex- 
amples of connections (for us) are the Levi-Civita connection of a Riemannian 
manifold, and a principal connection of a principal bundle. We will review these 
next. To do this we introduce the notion of an affine connection. 



Definition 2.4. An affine connection on a vector bundle tt : E — >• M is a mapping, 
V : X(M) © T(E) -»■ T(E) such that 
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(1) V /W (a) = fVxcr + gVya for f,g e C°°(M),X,Ye 3L{M),a G T(E). 

(2) V x ((7a + a 2 ) = V x <7 a + V x o 2 /or a a , a 2 G r(£). 

(3) V x (/<r) = X[f]a + /V x a /or / G C°°(M), X G £(M), a G r(£). 

Let M be a Riemannian manifold with a metric (•, -)m- A connection on TM is a 
subspace, H(TM) C TTM and the covariant derivative is a map : C(I, TM) — > 
TM. The covariant derivative uniquely defines an affine connection by the equiva- 
lence (Vx(^))(wio) — j£ where v(t) = Y(m(t)) and m(t) is the integral curve of X 
with the initial condition m |AM00[ §2.7]. There is a one-to-one correspondence 
between affine connections and covariant derivatives. 

On a Riemannian manifold there is a canonical affine connection which satisfies 
the following two extra properties: 

(1) V X Y - VyX = [X, Y] (torsion free). 

(2) (W Z X,Y) M + (X,W Z Y) M = Z[(X,Y) M ] (metric) 

for any X, Y, Z G X(M). We call this unique affine connection the Levi-Civita 
connection. 

Proposition 2.4. Let a G fl 1 (M). A torsion-free connection (such as the Levi- 
Civita connection) on TM induces the equivalence 

where we are viewing a as a linear function from TM to M. on the right hand side, 
and an element ofVL 1 (M) on the left hand side. 

Proof. Let X, Y G X(M). By the definition of |f viewed function on TM we 
find 



/ da 
\ dx 



(X),y\ (x) = (d TM a,h x(x) (Y(x))) 



where dxM is the exterior derivative on TM. Note that h x ^(Y(x)) is a vector in 
TTM over X(x) G TM obtained by parallel transporting the vector Y(x) over a 
curve with velocity X(x). In other words 



(a(x(t)),Y t ) 



t=o 

where x(t) is an arbitrary curve such that x(0) = x and ff| t _ = X(x). By the 
induced covariant derivative on curves of covectors we find 



da {x) y\ '-^ / Da W)) 



dx y " / v ' \ Dt 



Y i} ) + la(x), D)! 



t=o 



Dt 



t=o 



However, as the Y t is parallel along x(t) the second term is 0. Therefore 



9a (X),Y)(x^- /Da{m 



dx y n r \ Dt 



t=o 
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a -> Dt 



As Y Q = Y(x) we see that we can use the identity 4 {a, Y) = (j^, Y) + 1 • L>) 
to rewrite the above equations as 

(a(x(t)),Y(x(t))) - (a, D 



<>a (X),Y X 



dx / dt 



(Y(x(t))) 



As x G M is arbitrary, we can write this in terms of the affine connection V as 

' da 



i)x (X),Y)=C x (a,Y)-(a,V x Y). 
If we swap X and Y and subtract these equations we find 

0(X {X),y) - (^(Y),x) = C x (a,Y)-C Y (a,X) - (a,V Y X - V X Y). 



dx I \ dx 

By the torsion free property of the Levi-Civita connection, this implies 

^{X),y\ - (^(Y),x\ = C x (a,Y)-C Y (a,X) - (a, [X,Y]). 

However the right hand side is simply da(X, Y). □ 



Note that Proposition |2.4| presupposes that the exterior derivative is well de- 
fined. This is never a problem on finite-dimensional manifolds. However, infinite- 
dimensional manifolds require that we deal with some functional analytic concerns 
which we would like to skip |AMR09l see §7.4]. However, one can verify that the 
above theorem is true in the infinite-dimensional case if the exterior derivative 
exists. So we must add this existence assumption in the infinite dimensional case. 

Given a vector bundle tc : E — > M we may define the set of .E-valued forms 
Q 1 (M;S) := T(E) <S> f2 1 (M). This allows us to define the exterior differential of 
an E valued form by d(e (g) a) = e®da for an arbitrary e G T(E) and a G Q 1 (M). 
It is simple to see that using this definition for the differential of an ^-valued 
form proposition |2.4| holds for i?-valued forms as well. Additionally, the choice of 
a connection (perhaps not a Levi-Civita connection) on TM induces an intrinsic 
means of writing the Euler-Lagrange equations. 

Proposition 2.5. Let q(t) be a curve in a manifold Q and let L : TQ ->Kfc« 
Lagrangian. Then the following are equivalent 

(1) The curve, q{t), satisfies the Euler-Lagrange equations 

D fdL\ _ dL _ 



Dt \ dq J dq 

with respect to an arbitrarily chosen connection. 
(2) The curve, q(t), extremizes the action 

S= [ T L(q,q)dt 
Jo 

with respect to variations of q(t) with fixed end-points. 
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Proof. Before we begin we note that we have an equivalence of mixed partials: let 
q(t, e) be an embedding of a surface into Q (i.e. a deformation of a curve). Then 
we observe that 

Dq ( ( d 2 q \ \ D5q 



(7) 



v i ver 



\dtde J J Dt 



5S 



d_ 
de 



where q = || and Sq = ^. 

Let q : [0, T] <^-> Q be a curve in Q and let e) be a variation of q(t) with fixed 
endpoints. That is to say, q(t, 0) = and <?(0,e) = q(0),q(T, e) = g(T) for any 
e. Then we find 

»T 

L((q,q)(t))dt 

=o -'o 

,T q 

(dL, — (q(t),q(t)))dt 
'o <9e e =o 

" r ,<9L r x .dL D 

„ W ,i,) + < *'c; (,,)> * 

T ,<9L r x ,&L D 

T rT 



dL 



dv 



Sq 



dL 

t=0 Jo dq ' 
T .dL r . .DdL e 



5q)dt 



dq Dt dv 

As the variation is arbitrary, it must be the case that |p — (|^) = where 



□ 



3. KINEMATICS 

In this section, we return to the fluid-structure Lie groupoid G, which we dis- 
cussed briefly in Q. We introduce a vector bundle A, whose elements are pairs 
of body velocities and Eulerian fluid velocities, and we show that A is the Lie al- 
gebroid associated to the Lie groupoid G. Those who have not seen Lie groupoids 
before ought not be intimidated, as we will not require the reader to have heavy 
background knowledge on the topic. However, a knowledge of Lie groups will be 
helpful because many constructions on Lie groupoids parallel constructions on Lie 
groups. 

3.1. The fluid-structure groupoid G. Recall that the fluid-structure groupoid 
G is given by 



G : = {(&!, <p, b ) | b , &i e Emb(fi), <p G SDiff (f 



J b ) 



l)}> 
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where SDiff(R«^ , fsw&j) is the set of all volume-preserving diffeomorphisms from 
p^d bo into p^dbi- The groupoid G is equipped with a source projection 7r (&i, <p, b ) = 
b and a target projection 7r 1 (6 1 , <p, b ) = b\. Two element (b 2 , t/>, &i), (bi, (p, b ) G G 
such that the target of the second matches the source of the first element may be 
"composed" via the partial multiplication 

(&2,V>&i) ° { b i,(f,b ) := [b 2 ,ipo Lp, b ). 

We say that this is a "partial multiplication" because we can not neccesarily mul- 
tiply any element of G with another, but instead we require that targets and 
sources match in the appropriate way. Clearly this multiplication is associative 
and elements of the form [b, Id, b) serves as identity elements. Given these iden- 
tity elements, the inverse of (b\,(p,bo) G G is given by (b , I P~ 1 , b\). The associa- 
tive property, along with the existence of identities and inverses make G into a 
groupoid. Moreover, G is a smooth manifold, and the multiplication and inverse 
operations are smooth so that G is a Lie groupoid. 

An element (p, bo) G G can express the state of a body immersed in a fluid 
after a finite time, but how do we express infinitesimal motions? It turns out that 
the infinitesimal analogue of Lie groupoid G is the Lie algebroid A (see Figure [2]), 
which we now describe. 

3.2. The Lie algebroid A. To describe the configuration of the fluid-structure 
system, we begin by choosing an embedding bo G Emb(£>) which represents the 
location of the body at a chosen, initial time. The state of the fluid together 
with the structure can then be described by the element (bo, Id, bo) G G, where 
Id is the identity. As the fluid and structure are set in motion, they describe a 
trajectory g(t) = (b(t), <p t , b ) in G, where <p t is a diffeomorphism from p^d bo to 
fWb(t)- Note that the source of g(t) remains bo for all time, so that we obtain a 
curve in the source fiber over bo. The infinitesimal motion of the fluid-structure 
system is obtained by differentiating this curve at t = 0, and so we propose the 
following space as the space of infinitesimal motions: 

A={§ eTG | g(0) = Id b ,no(g(t)) = 

Interpreted physically we may denote an arbitrary element of A b by (b, b, u) where 
(b,b) G T b Emb(B) and u G X vo \(p^d b ) is an Eulerian fluid velocity which does not 
penetrate the moving body. In particular, because the "fluid part" of G consists of 
volume preserving diffeomorphisms which move the fluid, the vector field u must 
be divergence free and satisfies the no-penetration boundary condition, i.e 

(9) V-m = 0, and (b{b~ l (x)) - u{x),n(x)) = for all x G b(dw bo ). 

We now let A be the disjoint union 

A= |J A b 

feeEmb(B) 
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which is a vector bundle over r : A — > Emb(i3) called the Lie algebroid of G. 
We recall that we arrived at A by taking infinitesimal variations of the identity 
element (b, Id, b) G G. From this point of view, the relation between the groupoid 
G and the algebroid A is similar to the relation between a Lie group and its Lie 
algebra; see Figure [2] In fact this analogy extends furture because A is equipped 
with a Lie bracket on sections. A section of r : A — > Emb(i3) is a pair of maps 
(X B ,Xp^) : Emb(i3) — > A where X^y G 3£(Emb(i3)) and 3£r~ is a map which 
assigns a Eulerian fluid velocity field to each b G Emb(B). The bracket on sections 
of A is given by 

[(X B ,X^),(Y B ,Y^)](b) = ([X B ,Y B ](b),[X^(b),Y^(b)]). 

Lastly, there is one more structure on A known as the anchor map, 

p : (b,b,u) G A ^ (b,b) G TEmb(£>). 

This completes our description of A. Note that (as a set) A appears to be 
capable of expressing the fluid structure equations in the form of Q and ^ 
but may seem far removed from the realm of Lagrangian mechanics on tangent 
bundles. However, A will be shown to be the quotient of a tangent bundle of a 
configuration manifold by a Lie group action. This relationship with a tangent 
bundle of a configuration manifold is what we will use to connect A to Lagrangian 
mechanics. 

3.3. The bundle TQ(bo) / G(bo). Given an initial embedding bo G Emb(£>) we may 
set the initial condition of the fluid-structure system to be the identity element 
(b ,Id,b ) G G. As stated in the previous section, the physical evolution of the 
system is given as a trajectory g(t) = (b(t), cpt, b ) G G where the source of g{t) is 
60 for all time. In other words, if the initial configuration of the body is bo then 
the system evolves on the fiber 

(10) Q(b ):=n \b ) = {(b,<p,b )eG} 

(11) ={{bM I ^GSDiff(^ b0 ,R^)}. 

We may alternatively write elements of Q(bo) as pairs (b, <p) when we want to 
reduce clutter, or as triples (6, y?,&o) G G when we want to emphasize that Q(bo) 
is a subset of G. As trajectories on Q(bo) are all physically meaningful we may 
use Q(b ) as a configuration manifold for a fluid-structure system and we could 
implement Lagrangian mechanics on TQ(b Q ) in the standard way. However, we 
will postpone Lagrangian mechanics for now so that we may first describe the 
particle relabeling group for the fluid. In particular consider the Lie group 

G(bo) ■= n \b ) H vrr^o) = SDiff^J. 

We see that given an element (6, if) ~ (6, <p, bo) G Q(b ) and an element ip ~ 
(b ,ip, b ) G G(b ) we may use the groupoid multiplication to get the element 

(6, (f)-ip= (6, (p, b ) o (b , if), b ) = (b, (p o fa b ) = (b, yj o ip). 
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In this way, we obtain a right G(&o)-action on Q(bo) which relabels the particles 
of the fluid. Moreover, this action may be lifted to the tangent bundle TQ(bo) by 
taking tangent mappings of the action on Q(b ). In other words, for any ip G G(b ) 
we may define its action on TQ(b ) by the mapping 

(b, b, ip, ip) G TQ(b ) ^ (b, b, ip o ip, ip o ip) g TQ(b ). 

This group action on TQ(b ) allows us to define an equivalence relation where two 
velocites (b,b,(p,<p) G TQ(b ) and (b, V , ip', ip') G TQ(b ) are said to be equivalent 
if b = b' and there exists a ip G G(bo) such that ip 1 — (p o ip and tp' — tp o ip. We 
denote the equivalence class of a velocity (b,b,ip,ip) G TQ(b ) by [(6,6,(^,0)] and 
we denote the set of these equivalence classes by TQ(b )/G(bo). We won't bother 
with the technicalities to much, but suffice it to say that the set TQ(bo)/G(b ) 
is a manifold in its own right (see [AMOO, Prop 4.1.23]). Moreover the following 
proposition shows that TQ(b )/G(bo) is equivalent to the Lie algebroid A. 

Proposition 3.1. The map which sends Lagrangian to Eulerian velocities, given 
by 

^ : (6, b, ip, <p) i — y (6, b, u), where u = ip o ip , 

induces a Lie algebroid isomorphism between the quotient bundle TQ(b )/G(b ) 
and the Lie algebroid A. 

Proof. The space TQ(bo)/G(b ) is itself an algebroid, where the anchor map is 
given by [(6, b, <p, Lp)] !->■ (6, b) and the vector bundle projection is [(6, b, ip, ip)] i— > b. 
However, we will only show that ^ induces a vector bundle isomorphism; the fact 
that it sends the bracket on the space of section of TQ(b )/G(bo) to that on the 
space of sections of A is easily verified. 

The map ^ : TQ(b ) — > A given in the statement of the proposition is G(bo)- 
invariant, since 

q?(b, b, ip o ijj, (p o ip) = (b, b, ip o ip- 1 ) = q?(b, b, ip, ip) 

and hence drops to a quotient map from TQ(bo)/G(b ) to A, which we also denote 
by This map is linear on the fibers and is readily shown to be injective. All 
that remains is therefore to show that \1/ has an inverse which is also fiberwise 
linear. To construct this inverse, let {b, b,u) G A and choose an arbitrary ip such 
that (6, ip) G Qipo). Consider now the map $ which takes (b,b,u) G A into the 
equivalence class [(b, b, ip, u o ip)] G TQ(bo)/G(b ). A small calculation shows that 
$ does not depend on the choice of embedding ip, and that $ o = ldTQ(b )/G(b ) 
and $ o $ = Id^4, so that $ is the inverse of □ 



Remark 3.1. The isomorphism in (3.1 ) between the Lie algebroids A and TQ(bo)/G(b ) 
can be integrated to a groupoid isomorphism between the groupoid G defined pre- 
viously and the action groupoid (Q(b ) x Q(b Q ))/G(b ). 
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4. MECHANICS 

In this section we begin the physical description of fluid-structure interactions. 
The Lagrangian for fluid-structure interactions is the sum of two parts, one of 
which is fairly general and describes the immersed body, while the other is the 
kinetic energy of the fluid. Below, we will first explain the Lagrangian mechanics 
of a body in a vacuum before describing the mechanics of a body in a fluid. 

4.1. Elastic bodies in a vacuum. The Lagrangian for an elastic body is a func- 
tion Lb : TEmb(i5) — > R. We recall that the tangent bundle TEmb(i3) consists 
of all pairs (b, b) where b G Emb(i3) and b is a vector field along b; i.e. a map from 
B to TIR d such that the diagram 

TR d 
b ^ r „ 

B * R d 

b 

commutes. 

The kinetic energy is the function on TEmb(£>) given by 

K B (b,b) = l [ m(x)\\b(x)\\ 2 d 3 x 

where m : B — > R + is the mass-density of the body in the reference configuration. 
For an elastic body we could construct the potential energy by defining the right 
Cauchy-Green strain tensor, Cf, := &*(•, -) B . By construction is symmetric with 
respect to the metric on M. d and so we define the potential energy, Ui e (Cb), for 
some real valued function Ui e . Of course, any function U{b) on the shape of the 
body may serve as a candidate for the potential energy and choices of potential 
energy which differ from U\ e lead to various nonlinear theories of elasticity. In any 
case the Lagrangian for elasticity is 

(12) L B (b,b) = K B {b,b)-U{b), 

and the Euler-Lagrange equations of L B are the equations of motion for a body 
moving in a vacuum (see |MH83j for details). 

Remark 4.1. More generally, we may allow the Lagrangian of the body to be that 
of a gas, a rigid mechanical system, or even another fluid. We may even desire to 
add a surface energy, which would be particularly important if we were modelling 
bubbles or thin films (however we may opt to use Emb(<9£>) in the case of modeling 
a bubble). Again, the specifics of the Lagrangian for the body will not be addressed 
in this paper. 
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4.2. Bodies in fluids. If the body has configuration 60 G Emb(£>) at time t — 
and bi G Emb(i5) at time t — 1, then the configuration of the fluid can be described 
by a volume preserving diffeomorphism from R?y& to fsy^. This tells us that the 
configuration manifold for the fluid-body system with the initial body configuration 



b is precisely the manifold Q(b ) described in equation (11). We can define the 
kinetic energy of the fluid, : TQ(b ) — > K, by 

(13) W,^^) = - / IWaOfdvol. 

so that the total Lagrangian is given by 

(14) L{b, b, <p, ip) = L B {b, b) + Lp~(b, b, <p, 0). 

We may now "turn the crank" of Lagrangian mechanics to get equations of motion 
on TQ(bo). However, we will first reduce the description of the system by the 
particle relabeling symmetry to obtain equations of motion on A = TQ(bo)/G(bo). 

4.3. Particle relabeling symmetry. As in |Arn66] . the fluid-structure system 
is invariant under the action of G(b ) on Q(bo) by particle relabelings. 



Proposition 4.1. The Lagrangian, L, of equation (14) is symmetric with respect 
to the right action of G(bo) on TQ{bo). 

Proof. We observe that 

Lpaaib, b, cp o ip, tj> o ip) = I \\cp o -0|| 2 dvol 

bo 

\\<p\\ 2 detty" 1 ) dvol 



1 1 V? 1 1 2 g? vol 

where we have used both the change of coordinates formula and the volume preser- 
vation property, det(^) = 1. Therefore is invariant with respect to G(bo). The 
Lagrangian, Lg is symmetric with respect to G(po) in a trivial sense. Therefore 
the sum L = + Lg is G(b ) symmetric. □ 

For every symmetry of a Lagrangian system there is a conserved quantity via 
Noether's theorem. As in the case of classical hydrodynamics, we will find that 
the circulation of the fluid is this conserved quantity. 

Proposition 4.2. Let 0(60) denote the set of divergence free vector fields on &Vb 
which are tangential to the boundary. For each rj G 0(60) the circulation 

J v (b,b,(p,(f) = / ((p*r],(p o ip~ l )dvo\ 
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is conserved under the Euler- Lagrange equations of the Lagrangian, L, given in 
equation (14). 

Proof. A solution curve, q(t) = (b,(p)(t) G Q(b ), of the Euler-Lagrange equations 
extremizes the action 

S= L(b,b,ip,ip)dt 



with respect to variations with fixed end-points. Consider the variation 5q = 
S(b,ip) = (0,y? • rj) for r\ G $(b ). Since L is invariant with respect to G(b ) and 
we see that this variation comes from the deformation (b, <p o exp(rje)) we see that 
8S — 0. Additionally, we calculate that 

f T ,dL . . ,dL Dq\ , 

58 = I ^ 5 « ) + { ird )iL 

By the equivalence of mixed partials ^ this implies, 

f T ,dL . , ,dL D(6q), , 

T ,dL D fdL\ . ,dL r . " 



The first term vanishes because g(£) is assumed to satisfy the Euler-Lagrange 
equations. Using that 8S = then implies that (^,Sq) is a conserved quantity. 
This is Noether's theorem. We need only unpack the terms. We see that 



((p o (f 1 ,(p i ,rj)dvo\ 

b 

J v (b,b,<p,<p) □ 



Corollary 4.1 (Kelvin's Circulation Theorem). Given a closed curve C(t) G 
which is advected by the fluid and a vector field v along C(0), the quantity 

J v (b,b,(p,(p) — (f> (cp*v, (p o (p' 1 )d 3 x 



c 



is conserved by the Euler Lagrange equations of L. 



Sketch of the proof. This follows from Proposition |4.2| by choosing a sequence of 
vector fields which converge to one concentrated on C. See [AK92, Chapter 1] for 
details. □ 
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Given the G(&o)-symmetry of L we can see that the flow of the system exhibits 
a symmetry as well (as the flow is derived from L). To be precise, a curve q(t) G 
Qipo) which satisfies the Euler-Lagrange equations extremizes the action 

S[q] := f L(q,q)(t)dt. 
Jo 

For any ip G G(b ) we know that L((q,q) ■ if)) = L(q,q). This implies that the 
curve q(t) ■ ip G Q(b ) also extremizes the action so that q(t) ■ ip satisfies the Euler- 
Langrange equations as well. Therefore given one solution we obtain a collection 
of additional solutions; one for each way of relabeling the fluid particles via some 
ip G G(po). However, we are usually not concerned with the labels as they are 
chosen arbitrarily. In fact, because of the symmetry under particle relabelings 
there is a well defined flow on A = TQ(b )/G(bo) which is the space of elements 
(b,b,u), where the Eulerian velocity u is given u = (p o so that u does not 
depend on the choice of fluid labels. 

Because the Euler-Lagrange equations are derived from a Lagrangian, the G(bo) 
symmetry of the Lagrangian L should yield Euler-Lagrange equations which are 
also symmetric with respect to G(bo). Throughout the remainder of the paper we 
will move towards a description which expresses this idea and to do this we should 
recall how we have found the equations of motion in the past. In particular, it 
is standard to use Hamilton's variational principle to derive the Euler-Lagrange 
equations on TQ(bo)- This variational principal involved taking variations of curves 
in Qipo) and lifting them to variations in TQ(bo) by time-differentiation. We desire 
a variational principle on A = TQ(bo)/G(bo) which is consistent with Hamilton's 
principle, and so we must understand how variations of curves in Q(b ) are lifted 
to yield variations in A in an analogous manner. 

5. VARIATIONS IN A INDUCED BY VARIATIONS IN TQ(b ) 

We may derive the equations of motion on A in the same way that we derive 
equations of motion in Lagrangian mechanics, i.e. by extremizing the action in- 
tegral. However, taking variations of curves in A is not so simple. Given a curve 
(b, b, u) (t) G A a variation intuitively has two components: one component varies 
b, and the other component varies u. The issue with this perspective is that b 
is coupled to u through the no-penetration boundary condition Moreover, 
not all variations of u are physically meaningful, as we are only concerned with 
variations of u which come from variations of fluid particle motions (p. 

In this section, we will define the admissible variations on A as those which sat- 
isfy ^ and correspond to variations of fluid particle configurations. The language 
to to this will be aided by an isomorphism which decomposes these variations into 
"fluid" and "body" components. 
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5.1. Decomposing the variations in A. Hamilton's principle states that the 
action, J L(q, q)dt, is extremized along integral curves of the Euler-Lagrange equa- 
tions with respect to variations of curves with fixed end points. But what is a 
variation of a curve in this context? To answer this we first define a deforma- 
tion of a curve. Let I be an open interval on the real line containing 0. A curve 
q : t G / !->■ q(t) G Q is merely a mapping of I into Q. A deformation of the 
curve, q(t), is a mapping (t, A) G / x / h-» q(t, \) <E Q such that g(t, 0) = q(t). The 
variation with respect to this deformation is Sq := §f | A _ , which is a vector field 
along the curve q(t). The action is extremized with respect to arbitrary variations 
of q(t), but this does not imply extremization with respect to arbitrary variations 
of q :— ^. Specifically, the variations of q which we care about are those which 
are induced by the variation of q. In other words, we only care about variations 
of q of the form j- t 5q (i.e. we assume that q(t, A) is a smooth deformation). 

We have found that our system has a symmetry by the Lie group G(bo), and 
it is reasonable to assert the existence of a flow on the quotient space A = 
TQ(b )/G(b ). However, before we can explore this we must first understand 
the nature of variations of curves in Q(b ). Just as the variation of the curve 
(QiQ)(t) £ TQ(b ) is not arbitrary, we will find that admissible variations in A are 
not arbitrary either. 

Moreover, A can be seen, roughly speaking, as the product of two separate 
spaces, describing the body and fluid velocities, respectively. In order to express 
this splitting we may choose to decompose A into complementary vector bundles 
A = H(A) © V(A) where the vertical space is given canonically by 

V(A) :=kerne\(p) = {(b,0,0 & A}, 

with p : A — > TEmb(£>) the anchor map. In words, V(A) is the vector bundle of 
Eulerian fluid velocity fields which leave the body fixed. This vertical space comes 
equipped with the structure of a Lie-algebra on each fiber by the Lie bracket 
[(&, 0, £&), (b, 0, 77b)] := (b, 0, [£&, 77b]). We may define the adjoint map ad^ 6 : Vb(A) — > 
Vb(A) by ad^(6, 0,77b) = (6,0, [£&, 77&]). We call ad^ 6 the adjoint map, and we call 
the dual map ad^ 6 : V b (A)* — > V b (A)* the coadjoint map. The coadjoint map 
acts on the space of fluid momenta and will be important when we address the 
mechanical aspects of fluid-structure interaction. In summary, V(A) has a rich 
and canonically defined structure induced by A. In contrast, there does not exists 
a natural choice for the horizontal space H(A); we must choose one. 

5.2. The horizontal and vertical spaces. The only requirement on H(A) is 
that it is complementary to V(A) which means that we require p to be bijective 
on H(A). Therefore, one way to define an H(A) is by choosing a section X : 
TEmh(B) — > A which satisfies the property p ol = Identity. In other words, 
H(A) is defined by a section of p, denoted X, by setting H = range(X). Conversely, 
if we are given H C A, we may define the section X = (p\ H )~ l ■ 
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For X : TEmb(B) — > A to be a section of p means that for each (6,6) G 
TEmb(£>) there is a velocity field u given by (b,b,u) = 1(6,6). In particular, 
u G j£ vo i(f««b) is a velocity field which satisfies the no-penetration boundary con- 
dition Q since (6, 6, u) G A. Thus we may alternatively view X as a map which 
assigns (6,6) \-t u G X vo i(R^fe) such that m is a velocity field which satisfies the 
no-penetration boundary condition (see figure [3]). This identification of X as a map 
to Eulerian velocity fields will reduce clutter in future derivations and so this is 
how we will interpret X. 




Figure 3. The map X assigns to each body configuration 6 and 
body velocity 6 a velocity field u = X(6, 6) for the fluid, which is 
such that the no-penetration boundary condition holds. 



Remark 5.1. We note that there is considerable freedom in the choice of X, and 
as we shall see below, all choices of X yield different decomposition into fluid and 
body velocities but will still give the same equations of motion. The choice of X 
is analogous to the choice of a basis for an n-dimensional vector space. Different 
choices of basis yield different isomorphisms with IR n , and so a particular ODE 
on a vector space will be written differently on W 1 depending on which basis one 
chooses. 

We now make a small digression to show that one particular choice for X is 
singled out on physical grounds. For an inviscid fluid, we begin with the observa- 
tion that the vorticity is conserved throughout the evolution of the fluid-structure 
system. As a consequence, the Eulerian velocity field u for the fluid resulting from 
the body motion given by 6 is be a gradient flow, u = V$, where the potential $ 
satisfies the following Neumann problem: 

(15) A$ = 0, and —(x) = ({b o 6~ 1 )(a;), n(x)J for all x G b(dB); 

see |Lam45l IBatOO] . The solution of this equation is uniquely defined, up to a 
constant, and we now define a section of the anchor map denoted X by putting 



(16) 



X (6) = V$. 
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This section X can be viewed as the Eulerian version of the "Neumann con- 
nection" first introduced in [LMMR86J and further analyzed in |VKM0 9J. We will 
return to this observation in Appendix [AJ where we will show that there is in fact 
a one-to-one correspondence between the set of sections of p : A — > T Emb(B) and 
the set of principal connections in Q(bo). 

5.3. The covariant derivative on V(A). Let X be a section of the anchor map 
viewed as a map from TEmb(S) to Eulerian fluid velocities. Given a curve e h-> 
6(e) G Emb(£>) we may obtain an e-dependent sequence of vector fields 

u € =Z[b(e),- 

and ip £ will denote the flow of u e , i.e. ^P- = u e o ip e . We use the flow ip e to define 
a notion of parallel transport on V(A). 

Given a curve h e G Emb(i3) we may define the parallel transport of (bo, 0, £) G 
V(A) over b e to be the e-dependent curve (b e , 0, ((f e )*0 e V(A). This means of 
parallel transport allows us to define a covariant derivative. Using definition |2.1 



we see that one way to define a covariant derivative on the vector bundle V(A) is 
to define a horizontal sub-bundle to the vector bundle T(V(A)). We can do this 
by taking the derivative of the parallel transport to get an element 

(b e , 0, (if Me) = ((b, 0, 0, (Sb, 0, 5£ - [l(5b),Z])) 

=o 



(17 > * 



where Sb = ^| c _ n - The elements of the form (|17|) define a horizontal space on 



T(V(A)) and induce the horizontal and vertical projections 
hor((6, 0, 0, (6, 0, 0) = ((b, 0, 0, (6, 0, -[1(b), 
(18) ver((6, 0, £), (b, 0, £)) = ((b, 0, 0, (0, 0, £ + [l(b),£))). 



Proposition 5.1. The vertical projector given in (|18j) induces the covariant de- 
rivative 

(19) rJt {b ^ b)= ( b ' d a¥ + [I{h) ^ b 



Proof. This follows from Definition 2.2 for the covariant derivative. □ 



In Appendix |AJ we will show that this is in fact the formula for the covariant 
derivative induced by a principal connection (see |CMR0l|, IJRD12j ). 

Additionally the covariant derivative on V(A) induces a covariant derivative on 
the dual vector bundle to V(A) by equation ([6]). The dual bundle V(A)* is given 
by pairs (b, a^) where is a covector field on in the dual space to j£ vo i 
and the induced covariant derivative is given by 

(20) — — — = (b, a b - ad^ } a b ). 



24 



HENRY O. JACOBS AND JORIS VANKERSCHAVER 



5.4. The isomorphism between A and TEmb(B) © V(A). Now that we have 
thouroughly studied V(A), we can construct an isomorphism between the quotient 
bundle A and the direct sum TEmb(£>) © V(A). The advantage of this isomor- 
phism is that the summands in the latter space are easier to work with than the 
original quotient space. From a fluid-dynamical point of view, the isomorphism 
states the intuitively clear notion that the fluid velocity field can be decomposed 
into a part consistent with a moving body and a part consistent with a stationary 
body. 

Proposition 5.2. Given a section of the anchor map, X, the map ^% '■ A — > 
TEmb(S) © V(A) given by 

*i(6,6,«) = (b,b,u-X(b)) 
is a vector-bundle isomorphism which preserves the bracket structure on V(A). 

Proof. To check that ^% truly maps to the correct space it suffices to check that 
u — 1(6) is a vector field which is tangent to the boundary of pzvb- However this 
is true by construction because both 1(b) and u must satisfy the no-penetration 
condition given in equation ([9| and so the difference 1(b) — u is a vector field which 
is tangent to the boundary. Additionally, ^% is invertible with the inverse 

^(6,6,6) = (6,6,£ 6 +Z(6)). 

That \&z preserves the brackets on sections is proven by inspection. □ 



As an illustration, we write down the isomorphism in Proposition 5.2 in the 



case where the Neumann map defined in (16) is used in place of X. In this case, 
J (6) = V$. We then have that 

<po<p~ 1 = (<p o cp- 1 - V$) + (V$) 

consistent with stationary body consistent with moving body 



where $ is given by (15). The first factor, (fio(p~ — V$, is a divergence-free vector 
field tangent to the boundary, while the second factor is a gradient vector field, 
so that this decomposition is nothing but the Helmholtz-Hodge decomposition of 
vector fields. If we let P denote the projection onto divergence-free vector fields in 



the Helmholtz-Hodge decomposition, the isomorphism in Proposition 5^ becomes 

^([(6,6,^,0)]) = (6, 6, P^o^ 1 )). 

In general, when X is not a gradient vector field, the isomorphism \&x can be viewed 
as giving rise to a generalized Helmholtz-Hodge decomposition. 

5.5. Covariant variations. Let (b,ip) t be a curve in Q(b ). Then a deformation 
is a A-parametrized family of curves (6, (p)\,t such that (6, <f)o t t = (6, <p)t- We desire 
to measure how much the variation 5(b, cp) t = J| L=o ^' ^) A >* cnan g es the quantity 
(6, £b)t '■= (h, 0t°y?r 1— Zip)) To do this we invoke the covariant derivative 

induced by X of Proposition |5.1[ 
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Additionally, the kinetic energy metric, (-, -)Emb(B), on Emb(i3) induced by the 
kinetic energy, Kg, induces a Levi-Civita covariant derivative on Emb(£>). The 
Whitney sum of these covariant derivatives is a covariant derivative on the Whitney 
sum TEmb(B) © V(A) which we also denote 

Using this covariant derivative, we may define the covariant variation of a curve 
(b,£,b)t £ V(A) with respect to a deformation (b, £b)\,t by 

(21) ^,6)-^ 



D\ 



A=0 



We now compute the variation of a reduced curve (6, £&) f induced by a varia- 
tion of the curve (b,<p)t in Q(b ). We split this computation into two parts: in 



proposition 5.3 we compute the effect of vertical variations, which leave b fixed 



and act on ip by particle relabelings. In proposition 5.4 we then consider the effect 



of horizontal variations, which change b (and change ip accordingly). 

Proposition 5.3. Let I be a section of the anchor map of A and let (b,(p)t be a 
curve in Q(bo). Define the time- dependent vector field 



(€b)t = fit o <p t 1 - Z(b 



th 



so that (b, £b) t G V(A) for all t. Given a vertical deformation of of the curve given 
by (b t ,(pt°?pt,\) for a time- dependent deformation of the identity, ip tjX £ G(b ), the 
covariant variation of (6, £j,) t is given by 



5 X (b^b) = ^^-[(b^b),(b, Vb )] 



where 7] b = (tp t ) * V and rj := _ € g(b ) 



A=0 



Proof. To reduce clutter we shall suppress the t and A dependence of (p\,t and 

d_\ 

d\ I A=0 



abbreviate it as t ip\ The partial derivative 5£b '■= Jxk=n(£ fe ) ma y De S P^ 



three parts. 



^ dfi i m . dip' 1 dl(b) 

T x T 2 T 3 
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As the variation is vertical, 1(b) does not depend on A and therefore T3 = 0. We 
can rewrite T\ as 

O if 







T 1 = -(0)o^ 



dXdt 

|(M°^ 
— [v°v) 

Tup o 77 o y? -1 + Ty? o 77 o 
Tub or] o y? -1 + (^t), 



where we have used the fact that 

OA A=0 

On the other hand, T2 may be written as 

dip- 1 



Tip o 77. 



T 2 = Tip 



—T{ip o ip 1 ) o Sip o ip 1 

—T(ip o ip>~ 1 ) o Tip o 77 O <£> 
— o 77 o y? -1 



Thus we find 6£ b = y?*?). Since 56 = 0, and using the expression (21) for the 
covariant derivative, we see that 

5 x (b,£ b ) = (b,5£ b + [l(5b),£ b }) 



Additionally note that 



D(b, Vb ) 



Dt 



(6, (p*fj). 

6,^ + [1(6), 77,]). 



We calculate 



drj b d 
dt dt 



<p*r] = -[(p o (p 1 , ^77] + p>*r] 



Therefore ipj] = ^ + [f + 1(6), ^ so that 



^(6,£ 6 



+ X(6),77 fc ] + ^77. 

• th 

dr} b 



6,^ + [6+X(6),7/ b ] 



g(M>) 

Dt 



+ [(6,6), (6,^)]- 



□ 
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We understand how variations of curves in Q(bo) induced by G(bo) are expressed 
as variations in V(A), but how about other variations? We must also consider how 
{b,£b) — (b-><fi ° — X(fe)) varies in response to variations of b. Given a curve 
(b, ip)t G Q{bo) we may take a deformation of the curve bt G Emb(i3) given by bt,\ 
and define the horizontal deformation of (b, ip) t as follows. For fixed t, we consider 
the vector field u t \ defined by 



Ob 



d\ 



and we let (pt,\ be the flow of Ut,\, so that (pt,\ = Ut,\ o <pt,\- In this way, we obtain 
a family of curves (b tj \, <ft \), which we will call a horizontal deformation of (b, ip) t . 

Proposition 5.4. Given a horizontal deformation of(b,ip) given by (&t,A, ft,x) the 
covariant variation of (b, £&) := (b, up o ip~ x — 1(b)) is 

6 I (b,^ b ) = B(b,6b), 

where B is a V(A)-valued 2-form given by 



Bib, 6b 



dl b,6b)- 



1(b), X (6b) 



2.3 



and the fact that the fiber derivative, f|, is 

' ab ' 



Proof. Upon invoking Proposition 

identical to X (because X is linear in the fibers) we find that 

dip 



5& 



d\ 



-i 



A=0 



' dl 

Tip o ip' 1 o 6ip o ip^ 1 — ( ^ (&)> ^ 



X 



D 
Dt 



(6b) 



Using the equivalence of mixed partials we find that 



dip 
dX 



A=0 



d 
dt 
d 
dt 



{Sip) o ip 1 

(l(6b) o ip) o ip^ 1 



dj_ 
db 



(6b)^+l^ t (6b)^+l(6b)-ipoip~\ 



Upon substitution into the last line of the previous calculation we find 



6£ b 



m 

db 



(6b), b 



81 



db 



(b),6b ) + X(b, 6b) ■ (ip o v - 1 ) - (ip o ip" 1 ) ■ X(6b) 



and by proposition |2.4| 

= dl(b, 6b) + \X(6b),ubip- 1 } 
= dl(b, 6b) + [X(6b),i b + X(b)}. 
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Therefore the covariant variation is 

6 x {b,£ b ) = (b,8& - [Z{6b),&]) = (b,dX{b,Sb) - [z(&),Z(<y&)]) := B(b,8b). □ 

In summary, a variation Sq(t) of a curve q(t) G Q(bo) will lead to a variation of 
the curve q(t) G TQipo) which can be passed to the quotient A = TQ(bo)/G(b ). 
Moreover if we choose a section of the anchor map, Z, we get an isomorphism from 
A to the space TEmb(B) © V(A) and the relevant variations take the form 



Db Dr] b 



De Dt 



-[ Vb ,£ b }+B(b,6b) 



Now that we fully understand the journey from variations of curves in Q(bo) to 
variations of curves in T Emb(B)®V(A) we may use Hamiton's principle on TQ(bo) 
to derive the equations of motion on T Emb(B)®V(A). These equations are known 
as the Lagrange-Poincare equations. 



Remark 5.2. The U(.A)-valued 2-form B defined in the statement of Proposition 5.4 
can be interpreted as the curvature of Z when it is related to a principal connection. 
We will return to this point in Section A.3 , but for now we just refer to B as the 
curvature of Z. 



6. REDUCED EQUATIONS OF MOTION 

In this section we will derive the equations of motion on TEmb(i3) © V(A) 
using a reduced version of Hamilton's principle known as the Lagrange-Poincare 
variational principle. We will first derive the equations of motion in geometric 



form (Theorem 6.1) and we will then show that these equations are equivalent to 
the equations and ^ in coordinates (Theorem 6.2). 



6.1. The reduced Lagrangian. We recall from proposition 4T that the total 
fluid-structure Lagrangian, L, is invariant under the action of G(bo) on TQ(b ), so 
that we may define a reduced Lagrangian i : A — > K by 

i(b,b,u) = L(b,b,<p,<p), 

where (p and (p are such that u = up o <£> _1 . Explicitely, 

£{b,b,u) = L B (b,b) + \ [ \\u(x)\\ 2 dvo\. 

In the previous section we found that A is isomorphic to T Emb(£>) Q)V(A), given a 
section of the anchor map Z. In terms of Z the reduced Lagrangian on T Emb(B) © 
V(A) (which we will also denote by £) is given by 

(22) E(b,b,Z b )=L B (b,b) + ±J^ ||£ 6 +Z(6)|| 2 dvol,. 
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With this reduced Lagrangian we may transfer Hamilton's principle from TQ(bo) 
to the reduced space TEmb(B) © V(A) and derive symmetry-reduced equations 
of motion. 

6.2. The reduced equations of motion. Now that we understand how varia- 
tions of curves in Q(bo) induce variations in TEmb(i3) © V(A) and that we have 
an expression for the reduced Lagrangian £, we may transfer Hamilton's principle 
to the reduced space TEmb(B) © V(A). In the following theorem, which is a 
particular instance of the Lagrange-Poincare reduction theorem of [CMROl], we 
derive the equations of motion that result from this reduced variational principle. 

Theorem 6.1 (Lagrange-Poincare Theorem). Let X be a section of the anchor 
map p : A — > TEmb(£>) and let q(t) = (b,ip) t be a curve in Q(b ). Let (b, b, (,b)(t) 
be the induced curve in TEmb(i3) © V(A), where = <p o tp^ 1 — X{b). Finally, 



let L be the fluid- structure Lagrangian of equation (14) and let £ be the reduced 



Lagrangian of equation (22). Then the following are equivalent. 



'1) The curve q(t) extremizes the action 



S = L(q, q)dt 



o 



with respect to arbitrary variations of q(t) with fixed end points. 
(2) The curve q(t) satisfies the Euler-Lagrange equations 



Dt \dq J dq 

(3) The curve (b, b, extremizes the reduced action 

s= [ t(b,b,£ b )dt 



JO 

with respect to variations, 5b, of the curve b(t) G Emb(i5) with fixed end 
points and covariant variations of (b, £5) of the form 

S X (b, 6,) = + (b, Vb }) + B(b, 5b) 

for an arbitrary curve (b,7]b) t G V(A) above b t . 
(4) The curve (b, b, £&) t satisfies the Lagrange-Poincare equations: 

(23) ^A^)-f h =%B, 

( dl 



(24) irA^r)=-^h b )[ 





{-) 


Dt 


\db) 


° 


( dC 


Dt 





where is a two-form on Emb(i3) given by 5^(6,56) = (^,5(6,56) 
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Proof. The equivalence of (1) and (2) was established in 2.5 The equivalence of 
(1) and (3) follows fro m ta king a d eformation, (b, (p)t,\, of (b, (p) t with fixed end- 
points. By Proposition 5.3 and 5.4 we see that S x (b, = j^(b, r] b ) + (b, [£ 6 , 775] + 
B{b, 5b)) and 56 is a variation of with fixed endpoints. Therefore, if q(t) 
extremizes S with respect to arbitrary variations with fixed end points, then clearly 
(b, b, £b)t extremizes s with respect to the variations of the desired form. Conversely, 
variations 5b and 5 Z (6, £&) = -^(b,r]b) + (6, [?7&,£&] + B(b,5b)) produce arbitrary 
variations of (b, ip) t by the formula 5ip = (I(5b) + 77&) o {p. Thus (1) and (3) are 
equivalent. 

Now we will prove the equivalence of (3) and (4). By assuming (3) and taking 
variations we find 



= 5 



,b,£ b )dt 



\ ( —,—(b))dt 



— , -^z(5b) ) alt 



db Dt 



9i ,rn , A /9£ , 
T / d£ D \ I d£ D fd£^ 

W ~ ad ^ (M6 y + \Wb ~ m UJ + 56 } rft 

^ UJ " ad(6 ^ {wJ Ab ' Vb) / + \d-b-D-t{Fb) + ^ 56 } dt 



As 77b and 5b are arbitrary we arrive at (4). The above argument is reversible, and 
thus we have proven equivalence of (3) and (4) as well. □ 



Of course most people have opted to write the equations of motion for a body 
immersed in a fluid in the form of equations ([I]) and ([2]). While the Lagrange- 
Poincar'e equations appear different, they are merely same equations in different 
coordinates. As the proof of this equivalence is rather long, we have relegated it 
to Appendix |B"} 



Theorem 6.2. Under the assumptions of Theorem 6.1, the Lagrange-Poincare 
equations may be written as 



D (dL B \ dL B 
(25) ~D~t \~db~) db~ ~ 



du 

(26) - + u-V« = -Vp 

at 
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where u = £5 + 1(b) , Vp is a Lagrangian parameter which enforces incompress- 
ibility, and F p G T* Emb(i3) is a force on the boundary of the body given by 

(27) (F p ,5b)=f {p(x)n(x),{Sbob- 1 )(x))d 3 x. 

7. VISCOUS FLUIDS 

To add a viscousity to our formulation we need to alter our constructions in 
two ways. Firstly, we incorporate the no-slip boundary condition by using the 
subgroup oid 

G ns := {(&i, <p, b ) G G | if o b \ dB = h\ dB } 

where the source and target projections 7r , tti are simply the restriction of the orig- 
inal source and target to this subgroupoid. As before we choose a reference con- 
figuration b G Emb(i3) and define the configuration manifold Q ns (bo) = 7Tq 1 (6 )- 
The group of fluid particle relabeling symmetries is the group of diffeomorphisms 
of 

xttbo which keep the boundary point-wise fixed, 

G ns (b ) :=7T - 1 (6 )n7rr 1 (6o) 

= {■0 6 SDiff(?«y& ) 1 ip(x) = x,Wx G diWb }- 

The Lie algebra of G ns (b ) is the set 

0ns(M = {V £ X vo i(r^ 6() ) I T}(x) = 0, \/x G dptt bo } , 

and the Lie algebroid ^4 ns changes accordingly: ^4 ns consists of all triples (b, b, u) 
such that u(x) = b(b~ 1 (x)) for all x G dtw. Similarly, the bundle V(A ns ) consists 
of pairs (b, £) such that C, = on the boundary of the domain. It can be proven, 
as before, that *4 ns is isomorphic to TEmb(B) © V(A QS ). 

The second change that we need to make is the addition of a viscous friction 
force. For a viscous fluid with viscosity v > 0, the viscous force is a map F v : A — > 
A* defined by 

{F v (b,Vb,u), (b,Wb,v)} = v l trace(Vw T Vv)<ivol = — v \ (Au,v)dvo\. 

Below, we will also need the lift of F v to TQ(b ), denoted by F v and defined by 

(28) (Fvfav^ipjVp), (b,w b ,(p,w v )) = (F u (b,v b ,v v o ip' 1 ), (b,w b ,w v o y?" 1 )) , 

for all (6, v b , (p, v^), (b, w b , (p, w v ) G TQ(b ). 

Because of the isomorphism between A and TEmb(B)(BV(A), we can also view 
F v as a map from T Emb(i3) © V(A) to T* Emb(S) ®V(A)*, denoted by the same 
letter and given by 

(F v (b,v b ,£),(b,w b ,ri)) = v [ trace (V(X(u 6 ) + £) T V(Z(u> 6 ) + 77)) d vol. 
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Moreover, F u can be split as the direct sum of two parts, F v = F ns © f v , where 
F QS : TEmb(B) ®V(A) — > T* Emb(i3) is the boundary force which will enforce the 
no-slip condition and f v : TEmb(B) © V(A) —> V(A)* is the viscous force on the 
fluid. Finally, the force F v can be lifted to a force F u : TQ(b ) — > T*Q(bo) in the 
obvious way. These forces allows us to derive the equations of motion for a body 
in a Navier-Stokes fluid as shown in the theorem below. 



Theorem 7.1. Let L be the fluid- structure Lagrangian of equation (14) restricted 
to G ns and I : TEmb(B) © ^(*4 ns ) — > M the reduced Lagrangian written in terms 
of a section X : TEmb(B) — >■ A of p. Given a curve q{t) = (b(t),<p(t)) G Q(b ), 
let (b,b,£b)t be the induced curve in TEmb(B) © V(A ns ), i.e. such that £ b = 
ip o ip^ 1 — 1(b). Then the following are equivalent: 

(1) The curve, q(t), satisfies the forced Euler- Lagrange equations 



Dt \dq J dq 



where F v is the lift of the viscous force F u , defined in (28). 

The curve (6, b, £b)t satisfies the forced Lagrange-Poincare equations 

d_ fm\ _ m 

Dt \db) ~ db 



hB^ + F n 



d\{m) -" ad (^) {w b ) +fv 

where, 

(F m (v b ,£ b ),w b ) =u (n(x) ■ Vu(x),w b {b' 1 (x)))d 2 x 



(f»{vb,£,b),r}b) = -v j (A(X(uj) +£&),77&)dvol. 



(3) The curve (6, b, £ b ) t satisfies 

du 

— — \-u- \u = —Vp + u/au 
dt 

D_fdLs\_dLs 

Dt\ d h) db + 

where u — £& + 1{b), and F p is the pressure force on the boundary of the 



body given in Theorem 6.1 



Proof. The equivalence of (1) and (2) is merely a statement of the Lagrange- 
Poincare reduction theorem with the addition of a force. Assume (2) and we 
will prove equivalence with (3). We will deal with the vertical equation first by 
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understanding the force on the fluid. We find that the total viscous force is given 
by 

(F v (b, &,&), (b,w b ,T) b )) = / trace(Vu T V(?7 b + l(w b )))dvo\ 



—v / (Au, r/ + I(w b ))dvo\ 

+ v I (n(x) ■ Vu) • (X(w 6 ))rf 2 S 



Using the boundary condition , = w b ob 1 \ d r^ db , we may re- write the 

above force as 



(F v (v b , (w b , ri b )) =u (n(x) ■ Vu(x),w b (b 1 (x))d 2 S 

+ / (Au(x) i w b (b~ 1 (x)))dwo\ 



Following the same procedure as for the inviscid case (see "step 1" of Appendix 
B|) we arrive at the equation 



u b + ad;(w b ) = vAu. 

on R^fe. On IR d this becomes 

ou _ _ . 
— — h u ■ \u = — vp + vAu. 
at 

Additionally, if we denote the horizontal Lagrange-Poincare operator 

D fd\ d 



CVhoT - Dt\dbJ db 
we find that (just as in the inviscid case), £7\ r(^R?y) = i b B + F p , where 

(F p (b } b),w b } = / p(x)(n(x),w b (b~ 1 (x)))d 2 S 



C) r r 



is the pressure force on the body and 
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is the reduced Lagrangian for the fluid. Note that i = Lb + ipsa and CVhor 
is a linear operator on the set of reduced Lagrangians. Therefore the (forced) 
horizontal Lagrange-Poincare equation states 

Putting all these pieces together yields the forced Euler-Lagrange equation on 
Emb(i3) given by 

D_ (d_Ls\ _dLs_ 
Dt\db) db Pp + **" 
Reversing these steps proves the equivalence of (2) and (3). □ 

8. CONCLUSION 

In this paper we have articulated the Lagrangian picture of fluid-structure sys- 
tems as systems evolving on the source fibers of a groupoid, G. This was done 
in a way which was remarkably parallel to the observations of V. I. Arnold in 
regards to fluid systems as evolving on Lie groups |Arn66] . Moreover, we have 
shown that the equation can be right trivialized and thus there exist equations of 
motion on the Lie algebroid, A. Upon choosing a section X : TEmb(i3) A, we 
found that these equations are the standard ones obtained via the classical theory 
of continuum mechanics as in |BatOO] . We were then able to derive the equations 
of motion for a body immersed in a Navier-Stokes fluid by adding a viscous force 
and a no-slip condition to our groupoid. Given this geometric perspective we can 
extend, specialize, and apply our findings. Such endeavors entail the future work 
for this research. 

8.1. Rigid bodies. We could restrict our groupoid further by considering a sub- 
groupoid 

G ri gid := {9 e G | 3AeSE(3),f3(g) = A-a(g)}. 
We expect the resulting equations of motion to be that of a rigid body immersed 
in an ideal fluid. This would be a "cute" generalization of Arnold's discovery that 
the Euler equations are right trivialized geodesic equations on a Lie group, in that 
we would have found that the equations for a rigid body in an ideal fluid are 
right trivialized geodesic equations on a Lie groupoid. The equations of motion 
for a rigid body in a fluid have been treated from a geometric point of view in 
|Rad03l IKMRMH051 IVKM091 IVKMIO] , but so far the relevance of groupoids in 
this context has not been explored. 

8.2. Complex fluids. There exists a unifying framework for understanding com- 
plex fluids using the notion of advected parameters |GBR09j . This framework 
uses the interpretation of an incompressible fluid as an ODE on the set of volume 
preserving diffeomorphisms and extends the ideal fluid case by appending param- 
eters in a vector space which are advected by the diffeomorphisms of the fluid. 
This framework captures magneto-hydrodynamics, micro-stretch liquid crystals, 
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and Yang-Mills fluids. Merging these ideas with the groupoid interpretation of 
fluid structure interactions presented in this paper should be possible. Moreover, 
the notion of "advection" fits naturally in a groupoid theoretic framework where 
the source map stores the initial state of the system. 

8.3. Numerical algorithms. Discrete Lagrangian mechanics is a framework for 
the construction of numerical integrators for Lagrangian mechanical systems. These 
integrators typically preserve some of the underlying geometry, and as a conse- 
quence have good long-term conservation and stability properties (see |M WOlt 
HLW02J). The ideas of discrete mechanics have been extended to deal with me- 
chanical systems on Lie groups (see |LLM08] and the references therein) and in 
[GMP + ll] a variational integrator for ideal fluid dynamics was proposed by ap- 
proximating SDiff (M) with a finite dimensional Lie group. It would be interesting 
to pursue a melding of these integration schemes by "discretizing" the Lie groupoid 
Grigid to obtain a variational integrator for the interaction of an ideal fluid with a 
rigid body by approximating it as a finite dimensional Lie groupoid. The groupoid 
framework of this paper would provide a natural inroad into this problem, the more 
so since in |Wei95j it was shown that many variational integrators can viewed as 
discrete mechanical systems on Lie groupoids. 

Appendix A. PRINCIPAL CONNECTIONS 

The goal of this appendix is to show that the choice of a section 1 : T Emb(i3) — > 
A of the anchor map p : A — > TEmb(i3) defines a principal connection on the 
principal fiber bundle n : Q(b ) — > Q(b )/G(b ). The material described here 
relies on [KN63J; see also |CMR01j . This manifestation of a section of the anchor 
map appeared in early work on particles immersed in fluids where such an object 
was called an "interpolation method" [JRD12J. For the sake of simplicity, we treat 
the case of an inviscid fluid; the extension to viscous flows is straightforward. 

A.l. The connection one-form. Given a section 1 of the anchor map p : A — > 
TEmb(i3), we define a one-form A on Q(b ) with values in g(b ) by putting 

(29) A(b, (p) ■ (6, p) = Tip' 1 op-p* [1(b)]. 

In other words, A is a map from TQ(bo) to g(&o)- The right-hand side of this 
definition indeed gives rise to an element of g(&o); i-e. a divergence- free vector field 
which is tangent to the boundary of ^b - To see this, note that A can be written 

as 

(30) A(b, <p) ■ (b, p>) = <p* (> o p- 1 - 1(b)) , 

where the expression between parentheses on the right-hand side is an element of 
$(b). Since p is volume-preserving and maps the boundary of to that of 
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we have that the pullback vector field is an element of Q(bo). We may also write 
this correspondence as 



(3i) A(b, ^ ■ (&, cp) = ^ (e R (cp, if) - i(b) 

where 0# is the right Maurer-Cartan form, defined by Qr(p>, ip) = ipo if' 1 . While 
this expression is a bit more involved than the original definition, it will make the 
computation of the curvature of A below easier. Strictly speaking, Or is not a 
Maurer-Cartan form in the usual sense, but rather lives on the Lie groupoid G. It 
turns out, however, that Or has all the properties required of it — in particular, 
Or has zero curvature — so that we will not dwell on this point any further. For 
more information about groupoid Maurer-Cartan forms, see [FS08J. 



Proposition A.l. The one-form A defined in (29) is a connection one-form. 

Proof. We need to check two properties: 

(1) A is (right) equivariant under the action of G(bo) on Q(bo): for ip G G(b ) 
and (b, <p, b, ip) G Q(bo), we have 

A(b, cp o ip) ■ (b, cp o ip) = T(cp o ip)' 1 o cp o ip — (ip o ip)* [1(b)] 

= ip*(Tp>~ 1 o ip -1(b)) 

= Ad^-i [A(b,(p) ■ (b,<j>)], 

where Ad^,-i : Q(b ) — > Q(b ) is the adjoint action given by Ad,/,(£) = 
Tip o £ o 

(2) A maps the infinitesimal generator £q(& ) of an element £ G fl(&o) back to 
£. Recall that the infinitesimal generator is the vector field £q(6 ) on Q(bo) 
given by 

^Q(b )(b,<p) = ( Q > T( P°0 e T (b, v )Q(b ), 
for all £ G g(bo) to £, so that 

A(b, if) ■ £ Q( 6 ) = Tip- 1 o Tip o £ - 1(0) = £. □ 

Conversely, given a connection one-form A, we may define a section of the anchor 
map 1 as follows. For (b,b) G TEmb(i5), choose <£> and ip to be compatible with 
(b, b), that is, so that (b, <p\ b, ip) G TQ(b ). We then put 

(32) 1(b) = <pop- 1 - ( p* [A(b, ip) ■ (b, p)]. 

We first check that this prescription is well-defined, i.e. that 1(b) does not depend 
on the choice of p and p. Let ip' and dp' be different elements so that (b, ip'; b, p') G 
TQ(b ). There then exists a diffeomorphism ip G G(b ) and a vector field £ G fl(&o) 
so that ip' = p> o ip and ip' = <p o ip -\- T(ip o ip) o £. After some calculations, we then 
have that 

0' o i^' -1 = ip o y?" 1 + (<p o ■?/))„£ 
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and 

^[A(b, if') • (6, 0')] = <P*[AQ>, <p) ■ (b, 0)] + [ip o VU, 

so that 

^ ° V -1 - <P* [^4(6, V 2 ; 6, 0)] = V ° V 5 '" 1 - V 2 * [A{b, <p'\ b, <p')] , 

i.e. X is well denned. 

It is not hard to show that X satisfies the slip boundary condition ([9]), so that 
we arrive at the following result. 



Proposition A. 2. The map I defined in (32) is a section of the anchor map 
p: A^TEmb(B). 

A. 2. The induced horizontal bundle. The horizontal sub-bundle induced by 
a connection one-form, or equivalently, by a section of p is defined as follows. At 
every point (b,<p) G Q, we define a horizontal subspace H(b,<p) C T^^Qibo) by 

H(b,(p) = keiA(b, (p), 

and we let H(Q(b )) be the disjoint union of all of these spaces, so that H(Q(b )) 
forms a sub-bundle of TQ(bo). 

In our case, we have that the elements of the kernel of A(b, if) are the tangent 



vectors (b,(p) which satisfy by (32) that 1(b) = <fi o tp 1 . As a consequence, there 



exist horizontal and vertical projectors hor, ver : TQ(bo) — > TQ(bo) given by 

hor(6, if] b, (p) = (b, (p; b, 1(b) o ip) 
ver(6, ip; b, <p) = (b, <p;Q,<p- 1(b) o ip). 

The horizontal bundle is therefore spanned by vectors of the form (b, (p; b,X(b)o(p), 
where (b, ip) G Q and b G T^B are arbitrary. 

A. 3. The curvature of X. Using Cartan's structure formula, we have that the 
curvature of a (right) principal connection is given by the two- form B on Q(b$) 
with values in 0(60), given by 

B = AA - [A, A]. 



We recall from (31 ) that the connection is given by A = — X, and we use the 
fact that Qr has zero curvature: dO R — [Or, Or] = 0. From this, it easily follows 
that the curvature B can be expressed directly in terms of X by 

B = -dX+ [XXI. 



Up to a sign, this is precisely the curvature tensor that appeared in Proposi- 
tion (Kty. 
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A. 4. The induced covariant derivative. Lastly, we use the fact that a connec- 
tion one-form A induces a covariant derivative on the adjoint bundle q to justify the 



expression (19) given previously. It is well-known (see e.g. |CMR0l[ Lemma 2.3.4]) 
that the covariant derivative of a curve t i— > (&(£),£&(£)) G is given by 



D 
Dt 



dt 



-[A(b,(p) 



where (b, ip) G Q(b ) projects down onto b G Emb(£>), and (b,(p) G T^^Qipo) is 



Substituting the expression (30 ) for A in terms of the section 



such that ipoip 

X, we obtain the formula for the covariant derivative given in Proposition |5.1 



Appendix B. PROOF OF THEOREM O 



In this appendix, we provide a proof of Theorem 6.2 We assume that the 



Lagrange-Poincare equations (23) and (24) hold. Our goal is to show that these 



equations yield the same dynamics as equations (26) and (25). This proof is rather 



long and it is easy to get lost, so we will begin with a roadmap. 
Roadmap. 



STEP 1. We begin by proving that the vertical equations (24) are equivalent to the 



inviscid fluid equations (26). 



STEP 2. To prove the equivalence between the horizontal equations (23) and the 



body equations (25), we note that the horizontal operator 

= d_ (dtf} _ d(- 
dt 



CP 



hor I 



db 



Ob 



is linear on vector space of real valued functions of TEmb(i3) © V(A). 
Therefore the right hand side of the horizontal LP-equation may be written 
as 

(33) CV hoI (L B ) + CV hor (£^) = B^b, ■)• 

STEP 3. The main technical part of the proof consists of computing £"Phor(fe«)- 
To do this, we proceed as follows: 
step 3a. We first compute 

step 3b. We then compute ^ 



ah 



step 3c. We finally observe that the Lagrange-Poincare equations can be writ- 
ten as £7\ r(^R?y) = F + B^(b, ■) for a judiciously chosen F. 
STEP 4. By the linearity of CVhor, subtracting the previous computation from the 



Lagrange-Poincare equations (33) implies that 

CV hor (L B ) = F 

STEP 5. Lastly, we prove that F = F p , where F p is the boundary force defined in 



(27). This concludes the proof of the equivalence between (23) and (25). 
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Step 1. We begin with step (1) in this roadmap. It is simple to verify that 
dl 



<96 



■Mi 



(i b +i{b\5i b )dvo\-- 

20 we see that 



(u, 8£b)d vol. 



In other words, j^r = u . By equation 
If we substitute these identities into the vertical LP equation we get 



m + ad X(6)( M 



ad (6,6)( w )• 



which implies 



0. 



The last line is Arnold's description of the Euler equation ^+u- X7u = X7p |AK92| 
§1.5] and completes step (1) of the roadmap. 

Step 2. We can now move to step (2) which deals with the Lagrange-Poincare 
operator. By inspection we can see that £Phor is linear and so the horizontal LP 
equation can be written as 



CP{L B ) + CP{ir r 



BJb,-). 



where is the kinetic energy of the fluid, so that I = L B + £^~. Moreover, note 



that 



db db 



+ 



9Lb 
db ■ 



Step 3. This brings us to step (3a) of the roadmap, which is the computation of 



By definition ^-(b, b, &) e T* Emb(£) is given by 



96 



{b,b,Z),6b 



db 



de 



e=0 



where (b e ,b £ ,^b,t) is an e dependent curve in TEmb(i3) © V(A) such that 8b 



de I e=0 



b e and 



{be,be,£b,, 



0. 



By equation ( 19 ) this implies 

(34) 



d 
de 



:6,6 



-[J(<H0,&]. 



£ = 
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By the Reynolds transport theorem we find that 



e=0 



U btt +I(k)\\ 2 dvo\ 



(35) 



\\£ b +X(b)\\ 2 X{5b)-n b dA 
If d 



2 W de 



6=0 



\\Z b .+Z(b e )\\ 2 )dvo\ 



Here we have used the fact that the boundary of tw bt moves with velocity X(8b) 
as we vary e from 0. We can then focus on the second term in the above sum and 
apply the covariant derivative on T* Emb(i3) © V(A)* to find 



1 d 

2 de 



e=0 



U b , e (x)+X(b e 



Ux)+X{b){x) 



d 
de 



(^ e (x)+X(b e )(x)) 



e=0 



Z b {x)+X{b)(x),-[£ b ,X(6b)} + 



d 
de 



(4)) 



e=0 



&{x)+X{b)(x), -fa,I{6b)] + ( ^(b),Sb\ + X ( ^ 



where the second line is obtained from the first line via equation (34) and the third 



line is gained by viewing X as a vector-valued one-form and invoking proposition 



Because b e is a parallel translate of b we know that ^ = 0. This allows us to 



Dc 



drop one term in the last line. Substitution into equation (35) yields 



Ob 



■-(b,b,Z),Sb 



\\i h {x)+X{b){x)\\ 2 X{5b)-n b dA 



+ 



Ux)+l(b)(x), -[£ b ,X(5b)} + (—(b),5b))dvo\ 



3X-, 



Finally, we may reduce clutter by substituting u — + X(b) and invoking the 
divergence theorem to transform the surface integral into a volume integral. This 
yields 



(36) 



db 



ib,b,£b),8b 



div ( ^\\u\\ 2 X(5b) 



+ (u,[X(8b)^ b } + 



dX 
db 



d vol 
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This completes step (3a) and we may move to step (3b). We find that the fiber 
derivative of is given by 

^{b,b,ti b ),5b 



db y ' ,w " / de 



£x~(b,b + e5b,£ b ) 
= (u\X(8b)). 

By the definition of the covariant derivative on V(A)* we find 
(37) (°(^), a \ /St ™ DSb 



Dt\ db dt\ \ db \ db Dt 



Ti T 2 



We can see by direct computation that 

This is as far as we need to go with T 2 and we will now rework T\. Firstly, we 
note that T\ is the time derivative of an integral over a time dependent domain, 
so we must invoke the Reynolds transport theorem a second time. This yields the 
equivalence 

Tx = (Z b + I(b),I(5b))(£ b +I(b),h b )dA + j t {t b + l(b),I(8b))dvo\ 
By construction (£ b , h b ) = on the boundary of p$n b so that 

71=/ (£ b +l(b),I(5b))(l(b),n b )dA+ [ l(u,I{6b)))dvoL 

Finally substituting u = 1(b) + £ b yields 

T x = (u,l(5b))(l(b),n b )dA + ( — ,l(5b)) + (u,-l(5b))dvol 



2.3 



a second time. Ex- 



The term |X(<56) can be handled by invoking proposition 
plicitly, this give us the equivalence 

where we have used the fact that X is fiberwise linear on TEmb(£>), and therefore 
X is equal to its own fiber derivative. We can substitute the above computation 
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into Ti to get the final expression 



Ti = J (u,X(5b))(X(b),n b )dA 



We can finally revisit equation (37), and write it as 



(u,X(Sb))(X(b),h b )dA 

+ L(f' iw ) + ("'(f w ^» dvo1 



We can substitute the vertical equation % = Vp — it • Vw and (as is customary) we 
will convert the surface integral into a volume integral via the divergence theorem. 
This yields 



d rat 



Sb^ = {div [{u,X{5b))X{b) 



+ {Vp- it ■ Vu,X{8b)) + /it, (^( Sb )^)) } rfvo1 

This bring us to step (3c) where we calculate £Phor(^R^)- By subtracting (36) 
from ( 38 ) we find 

m (1?) " ^ = L, l div O^^) - 5W^»)) 

+ (Vp - m • Vu,Z(6b)) + (u, [1(56), 6] > 
+ (*&),&)-(§ (b),5b\}dvo\ 



We may now invoke proposition 2.4 to replace the term ^||(<5&), b j — ^fj|(&) 

with dX(b,8b). Finally we use the equivalence Bib, 5b) = dX(b,5b) + [X(b),X(8b)] 
to rewrite the last equation as 

D fd£~v\ d£^ 



| div f(u,X(8b))X(b) - ^\\u\\ 2 X(5b) J 



Dt \ db J db Jr^ b i V 2 7 

+ (Vp - « ■ Vu,X(8b)) + (u, [J(<f&),&]> 

+ (it, [J(<J6),J(6)]> + (it,S(6, 56)))dvol 
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We define the term B^(b,5b) = fr^ b (u, B(b, Sb))d vol, so that we may conclude 
(39) CV har (£ f ^) = F + B li (b,-). 

Step 4. This brings us to step (4). Given that CP^ is linear on the set of 
Lagrangians the horizontal LP equations give us 

B,(b,-)=CV hoT (£) 

= £Phorfe) + CV hoi (L B ) 

= B fi (b,-) + F + £V hoT (L B ), 

so that 

CV(L B ) = F. 

Step 5. This brings us to step (5) where we will manipulate F. To reduce clutter, 
set v = I (5b) and w = 1(b) so that the term given by F may be rewritten as 

1 



L div (' 



u • v)w — -(u • u)v ) + u ■ [v,u] + (Vp — u ■ V«) ■ v d x 



Using Einstein's summation notation we may write F as: 

F = I — — u'v'w 1 " u i u i v k + u k [ -^-v' - --^v? + — ^-v k - u i v k ' 



... j dx k \ 2 / \ dx l dx % J \ dx k dx % 

By applying the divergence free properties, |^ = |^ = 0, this simplifies to 



-I 

= / div (p(x)w(x) — (w(x) ■ u(x)) £b( x )) d 3 x 

(p(x)I(5b)(x) - (X(5b),u(x))i b (x),h(x))d 2 x. 



Noting that = u — 1(b) = u — v is tangent to the boundary, and therefore 
orthogonal to the unit normal n we can drop the term involving so that 

F= [ (p(x)I(5b)(x),n(x))d 2 x 
and by the no-slip condition of X we find 

(p(x)n(x), (5b o b~ 1 )(x))d 2 x, 



so that F = F p , as claimed. 

Thus we may conclude that the Euler-Lagrange equations for L B are given by 

{w t (If ) " °w •») = L ° »"><*» ^ 
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for arbitrary variations 5b. Since the variation 5b is arbitrary the theorem follows. 
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